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Abstract 

We describe two-flavor QCD lattice data for the pressure at nonzero temperature and vanishing 
chemical potential within a quasiparticle model. Relying only on thermodynamic consistency, the 
model is extended to nonzero chemical potential. The results agree with lattice calculations in the 
region of small chemical potential. 
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I. INTRODUCTION 



One of the fundamental issues which triggered, and has influenced since, heavy ion physics 
is the question of the phase structure and the thermodynamic properties of strongly inter- 
acting matter at energy densities above lGeV/fm 3 . Under such conditions, exceeding the 
energy density in nuclei but still far away from the asymptotic regime, the coupling strength 
a s is large, which makes the theoretical description of the many-body problem challenging. 

In the recent past the understanding of this field has become much more detailed. The 
phase diagram for QCD with rif = 2 massless flavors, which is the case we will consider 
in the following, can be briefly described as follows (we refer to for a recent review). 
At zero quark chemical potential, fi = 0, the broken chiral symmetry of hadron matter 
is restored within the quark-gluon plasma, at a critical temperature T c « 170 MeV. It 
is thought that this second order transition persists also for nonzero /i, thus defining a 
critical line, which changes to a first order transition line at the tricritical point. For small 
temperatures and fi ~/i c one anticipates a color-superconducting phase of quark matter. 
The value of [i c is expected to be 100. ..200 MeV larger than the quark chemical potential 
fi n = 307 MeV in nuclear matter. Quantitative results for large a s can be obtained from first 
principles by lattice calculations which were, however, restricted to nonzero temperature and 
H = until very recently. Therefore, the described picture for ^ is mainly based on 
general arguments combined with results from various models, including extrapolations of 
perturbative QCD. 

As a phenomenological description of the thermodynamics of deconfined strongly inter- 
acting matter we proposed a quasiparticle model 

BB 

. Its parameters are fixed by the lattice 
data at fi = 0. We then use the fact that within the model the thermodynamic potentials 
at zero chemical potential and /i ^ are related by thermodynamic consistency. In [3J] we 
analyzed lattice data for rif = 2 flavors [J, and rif = 4 which were, however, still dero- 
gated by sizable lattice artifacts which have an effect on the absolute scaling of the data. We 
therefore introduced a constant effective number of degrees of freedom of the quasiparticles 
as an additional model parameter to obtain first qualitative estimates. Later we considered 



in [f| the lattice data Q, where also the physical case of (2+1) flavors was simulated. As 
the absolute scaling of the lattice data enters as important information in particular near 
T c , we pragmatically applied the continuum extrapolation of the data, which was proposed 
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in j?| for T > 2 T c , also for smaller temperatures. The results of this prescription can now 
be compared to new lattice data ^]. Meanwhile, there are other lattice calculations which 
allow us to test directly the assumptions underlying the quasiparticle model as well as, for 
the first time, some of its predictions for nonzero chemical potential. 

We will therefore consider here the presently available lattice data for rif = 2. Based on 
that, we will fit and discuss the quasiparticle parameters at fi = in Section 3. In Section 
4, we will briefly summarize how to extend the model to nonzero chemical potential, and 
compare our findings with the results [9] from lattice simulations studying the region of 
small /i. Section 5 concludes with the discussion of some physical implications. 



II. FINITE TEMPERATURE LATTICE DATA 



The simulations jq| are performed on lattices with spatial extent N a = 16 and temporal 
sizes N T = 4 and N T = 6, with an improved Wilson quark action and renormalized quark 
masses corresponding to fixed ratios m ps /m v of the pseudoscalar to vector meson masses. 
We first consider the data for two light flavors, corresponding to 0.6 < m ps /m v < 0.75. 
Although this is larger than the physical value, the results are almost insensitive to the 
ratio, which suggests that they are not too far from the chiral limit. As expected for the 
rather small lattice sizes, the results for N T = 4 and 6 differ. However, we observe that 
normalizing the pressure data by pQ° nt /p Q T , the ratio of the free limits in the continuum and 
on the lattice, improves considerably the consistency between the data sets. As a matter of 
fact, the normalized N T = 4 data are in agreement with the normalized N T = 6 data after 
rescaling by a constant of 1.14. This simple scaling behavior for large coupling is rather 
remarkable. Based on this observation we suggest the continuum estimate for the pressure 
shown in Fig. ^ We assume here that the normalized N T = 6 data are already close to the 
continuum limit. This is supported by the fact that the thus interpreted data match the 
aforementioned continuum estimate from the staggered quark simulations j]]] 1 . Therefore, a 
consistent picture forms for the thermodynamics of QCD with nf = 2 light flavors. 

In Fig. |2l the corresponding data for the entropy are shown. It is noted that since the 
slope of the continuum extrapolated pressure Q is slightly larger than that from the data 

1 In these calculations m q — 0.1 T was assumed, corresponding to m ps /m v = 0.7 at T c . From the weak 
quark mass sensitivity observed in [8j, both results should indeed be comparable. 
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FIG. 1: Compilation of rif = 2 lattice data for the pressure in units of the free pressure vq. 
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Shown are the scaled (see text) data |8| for light quarks corresponding to meson mass ratios of 
0.65 < m ps /m v < 0.75 (open circles: N T = 4; open squares: N T = 6), and the continuum estimate 
[3] (grey band). The full line is the quasiparticle result. The full symbols represent the data ^| 
for !ar g e q oar k masses, with ^ - 0.95 For comparison, the hatched hand shows the 80(3) 
lattice data (dotted line: [10(; dashed line: [llj) normalized to the corresponding free pressure. 



see Fig. [IJ, the upper part of the error band is already for T ~ 3T C very close to the 
free limit. This would be in contrast to the pure gauge case, where the uncertainty due to 
lattice artifacts has become small, so we will assume that the lower side of this estimate is 
more relevant. 



III. QUASIPARTICLE MODEL 



For completeness, we briefly recall here the main ideas of the quasiparticle model 




of the QCD plasma. 

For weak coupling g, the thermodynamic behavior of the system is dominated by its 
excitations with momenta ~ T. While hard collective modes (the longitudinal plasmon and 
the quark hole excitation) are exponentially suppressed, the transverse gluons and the quark 
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FIG. 2: The lattice data for the entropy corresponding to the data for the pressure shown in Fig. ^ 
and the quasiparticle fit. 



particle excitations propagate predominantly on simple mass shells, u)f{k) 
In the chiral limit the so-called asymptotic masses are given by 

1 



mf + A; 2 
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(1) 



where fi q denotes the quark chemical potential, and N c = 3. Interpreting the relevant 
excitations as quasiparticles, the thermodynamic potential is 



(2) 



where Pi = ±d{T J d 3 k/(2n) 3 ln(l ±exp{ — {yji — /ij)/T}) are the contributions of the gluons 
(with vanishing chemical potential) and the quarks (for the antiquarks, the chemical poten- 
tial differs in the sign), and d g = 2(N% — 1) and d q = 2N C count the degrees of freedom. As 
shown in [ljj, thermodynamic consistency requires the derivative, with respect to the m|, 
of the right-hand side of Eq. (j2j) to vanish, i.e., the contribution B is related to the T and 
\i dependent masses by 



dB _ dpj(^\jMjjn^] 



dm 2 - 



(3) 



This implies that the entropy and the particle densities are simply given by the sum of the 
individual quasiparticle contributions, 



dpi(T, fii,mf 



dT 



dpi(T, Hi;mf 
m = - 



dfj-i 



(4) 



while the energy density has the form e = J2i e i + B. 

Expanded in the coupling g, the above approach reproduces the leading-order pertur- 
bative results. The full expressions, however, represent a thermodynamically consistent 
resummation of terms of all orders in g. This suggests pondering the application of the 
model also in the strong coupling regime. 2 Considering first the case \l — 0, it indeed turns 
out that the lattice data for the entropy shown in Fig. |2] can be described by the model with 
the ansatz 

a ' (T ' = 0) = (lliV c -2rv)ln[A(r-r a )/r c ]2 (5) 
for g 2 /(4ir). This is the leading order perturbative result at a momentum scale determined 
by the temperature: T c j A is related to the QCD scale A, while T s parameterizes the behavior 
in the infrared. For the parameters we obtain 3 

A = 17.1, T S = 0.89T C . (6) 

The resulting quasiparticle masses are large; near T c they reach several times the value of 
the temperature. 4 The existence of such heavy excitations, which we have inferred from the 
thermodynamic bulk properties, has meanwhile been confirmed directly by lattice calcula- 
tions of the propagators "v\ . Finally, since the derivative of the 'bag' function B is related 
to the quasiparticle masses by Eq. (j3J), the model is completely defined by fixing 

B = B(T C ) = 1.1T C \ (7) 

which enters the fit in Fig. ^ as the third parameter. 



2 A formal reason supporting this attempt is the stationarity of the thermodynamic potential with respect to 
variation of the self-energies around the physical value; see 

14] 

and the references given there. Moreover, 
there are heuristic arguments that resummation improved leading order results might be more appropriate 



at large coupling than high order perturbative resu 
3 For the fit we considered only the normalized data 





The result then reproduces the extrapolated data 
on the lower side of the estimated error band; see the remark at the end of the last section. 
4 In an alternative approach, instead of attributing the deviations from the free limit at smaller temperatures 
to the mass of the quasiparticles, a variable number of degrees of freedom is proposed in 16]. 
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Since all the information about the coupling is encoded in the parameters T s and A, it is 
interesting to look at their flavor dependence. Comparing to the pure gauge plasma, it is 
recalled that in this case the pressure becomes very small close to the transition since there 
it has to match the pressure of the heavy glue balls in the confined phase. Similarly, the 
entropy is small at T ~ T c , which implies a large coupling there. For n/ = 2 the scaled 
entropy for T ~ T c is somewhat larger, thus close to the transition the coupling has to be 
smaller than for pure SU(3). However, for fixed parameters A and T s , the coupling (jSJ) would 
increase with increasing number of active flavors. Therefore, a difference of the parameters 
for n/ = 2 to those for the pure gauge plasma 0], 

ASU(3) = 4.9 ; T |U(3) = q 73 ^ ^ ^ 

is not unexpected. Interestingly, the parameter T s does not change by much compared to 
the case of rif = 2. 



IV. NONZERO CHEMICAL POTENTIAL 



The quasiparticle model as applied in the previous section can be generalized to nonzero 
quark chemical potential fi q = fx. The quasiparticle masses now depend also on [i - explicitly 
by the dimensionful coefficients of the coupling in Eq. ((TJ), and implicitly by the coupling 
itself. As shown in j^j], Maxwell's relation, ds/dfi = dn/dT, directly implies a partial 
differential equation for a s (T,[i). It is of first order and linear in the derivatives of the 
coupling (but nonlinear in a s ), 

where the coefficients ct, c m and C depend on T, /i and a s . It can easily be solved by 
reduction to a system of coupled ordinary differential equations, 

dT(s) _ dfi(s) _ da s (s) _ 

~ls~- CT ' ~Js~~ C ^ ~1s~- C ' (10) 

which determines the so-called characteristic curves T(s), /i(s), and the evolution of a s along 
such a curve, given an initial value. 

With regard to the underlying physics it is worth pointing out some properties of the 
flow equation Q. The coefficients are combinations of products of a derivative of the 



quasiparticle entropy or density with respect to the quasiparticle mass, and a derivative of 
the quasiparticle mass with respect to T, /x or a s . Writing down the explicit expressions, it 
is easy to see that the flow equation is elliptic. In particular, one finds 



The coefficient c M , e.g., vanishes because not only the entropy goes to zero as T — > 0, but 
also its derivative with respect to the mass. Therefore, the characteristics are perpendicular 
to both the T and the /x axes. This guarantees that specifying the coupling on some interval 
on the T axis sets up a valid initial condition problem. From the temperature dependence 
of the effective coupling as obtained from the lattice data at /x = 0, e. g. in the physically 
motivated parameterization (|5ji . we can therefore determine numerically the coupling from 
Eq. (J9|L and hence the equation of state, in other parts of the /xT plane. 

It is instructive to consider the asymptotic limit, a s — > 0, of Eq. Q, where the coefficient 
C vanishes. Then the coupling is constant along the characteristics, which become ellipses 
in the variables T 2 and /x 2 , leading to the mapping 



This holds approximately also for larger coupling, see Fig. |31 so the lattice data at /x = are 
mapped in elliptic strips into the /x T plane. On the other hand, an ansatz analog to Eq. (J5J) 
to parameterize a s (T = 0,/x) is quantitatively less satisfactory than in the case /x = 0. A 
closer look at the characteristics emanating from the interval [T c , 1.06T C ] reveals that they 
intersect in a narrow half-crescent region, which indicates that there the solution of the flow 
equation is not unique. This, however, is only an ostensible ambiguity. It so happens that 
the extrapolation of the pressure becomes negative in a larger region, see Figs. H3 and 0] This 
implies that a transition to another phase, at a certain positive pressure, happens already 
outside this region, so the encountered ambiguity of the flow equation is of no physical 
relevance. 5 

At this point we emphasize again that this extrapolation of the quasiparticle model relies 
only on the requirement of thermodynamic consistency. Of course, it implicitly assumes also 

5 We remark that the region where the solution of the flow equation is not unique is determined only by 
a s (p — 0, T), i. e. by the parameters A and T s fitted from the entropy, whereas the p = line depends also 
on p(p = 0,T C ) and thus on the third parameter Bq. Therefore, the fact that the potential ambiguity is 
irrelevant is based in a nontrivial way on the underlying lattice data for the equation of state at fx = 0. 



C T (T,/x = 0) = 



c M (T = 0,/x) = 0. 



(11) 




(12) 
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FIG. 3: Represented by the full lines are the characteristics of the flow equation 0. The char- 
acteristic through T c coincides for small u with the critical line (dashed, with a hatched error 
band) obtained in the lattice calculation 9]. In the region under the dash-dotted line the result- 
ing quasiparticle pressure is negative - a transition to another phase has to happen somewhere 
outside. Therefore, the narrow grey region under the p = line, where the solution of the flow 
equation is not unique, is physically irrelevant. Indicated by the symbol (assuming, for the scaling, 
T c = 170 MeV) is the chemical potential fi n in nuclear matter. 

that the quasiparticle structure does not change, i. e., that deconfined quarks and gluons are 
the relevant degrees of freedom. For small enough /i and temperatures above (or near, as 
(j, gets larger) T c this is a justified assumption. However, the quasiparticle structure will 
change in the hadronic phase, when both T and \i are small, as well as for sufficiently 
cold and dense systems where the color-superconducting phase is expected. Although the 
present quasiparticle model cannot make any statements about these phases, it is interesting 
to observe that it 'anticipates' the existence of another phase only from the lattice input 
at T > T c and \x = 0. An interpretation of the apparent similarity of the line of vanishing 
pressure in Fig. El with the expected transition line from the hadron to the superconducting 
quark matter phase, see Ref. jl|, remains, of course, a speculation. 

There is, however, a related question which we can address with the quasiparticle model 

9 



FIG. 4: The pressure scaled by the free pressure Po(T, /i); T and [i are in units of T c . The pressure 
along the characteristics starting out from T ~ T c becomes negative at small T, see also Fig. |31 
The change to a different phase has to happen already outside this region. 



without knowing details about the other phases, just based on the fact that for nonzero 
chemical potential the transition from the deconfined to the confined phase occurs at the 
critical line T c (/x). The critical line is expected to be perpendicular to the T axis, which has 
been confirmed in a recent lattice calculation :q| where also its curvature at /J, — has been 
calculated 6 , T c d 2 T c ({i)/d{i 2 \^ =0 ~ —0.14. Within the quasiparticle model it is natural to 
relate, at least for small /i, the critical line to the characteristic through T c (fi = 0), which, as 
shown above, is also perpendicular to the T axis. For small \x where only the quadratic terms 
are relevant (practically even for /x as large as 2 T c ) , we indeed find the T c characteristic in a 
striking agreement with the critical line from [j| , see Fig. 01 Another argument supporting 
the above interpretation of the T c characteristic comes from considering the case where the 
quark flavors have opposite chemical potentials, \l u = —[id — P>- With this isovector chemical 
potential the fermion determinant is positive definite, and standard Monte Carlo techniques 



181 ] - The lattice result H obtained for 



can be applied to study this system on the lattice 
the curvature of the critical line in that case agrees with the value quoted above for the 
isoscalar potential //. Within the quasiparticle model, the equality of these two numbers is 



6 In passing we note the amusing fact that the result agrees with the value from the bag model assuming 
free massless pions for the hadronic phase. 
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immediately evident . 

In Ref. |9| it was furthermore mentioned that the quadratic behavior, with the same 
curvature as at = 0, of the critical line is not likely to extrapolate down to small transition 
temperatures since T c (/z) would then vanish only at fi c ~ 650 MeV. Phenomenologically, 
however, fi c is expected to be not very much larger, say at most by 200 MeV, than the 
quark chemical potential fi n = 307 MeV in nuclear matter. In the quasiparticle model, from 
the chemical potential where the extrapolated pressure vanishes at T = 0, we estimate 
/i c «3T c ~ 500 MeV. 

This value is in the expected ball park, which encourages us to consider the extrapolation 
of the model down to smaller temperatures. Although for T — > quark matter will be in 
the superconducting phase, it is still possible to give an estimate of its equation of state in 
that region from the quasiparticle model. The quark pairing influences thermodynamic bulk 
properties at the order of (A fi) 2 , with the gap energy A being at most 100 MeV Q]. This 
has little effect on the energy density e = J2i e i + B as both the quasiparticle contributions 
and the function B are parametrically of the order 0(fi 4 ). For the pressure, on the other 
hand, the pairing effects become comparable to our expression p = J2iPi — B only when the 
latter becomes small. Since the pressure of the thermodynamically favored superconducting 
phase is larger than that of the plasma phase, the relation e(p) as shown in Fig.EJis therefore 
an upper estimate of the equation of state of cold quark matter. For p > 5 T 4 , we obtain 
e(p) m 13 T c 4 + 3.2 p, where the slope is mainly determined by the fact that the pressure at 
T = essentially scales as /i 4 . For smaller pressure, the slope is only slightly larger, and the 
energy at p = is approximately 7 11 T 4 . Assuming T c « 170 MeV, this translates into an 
energy density of 1 GeV/fm 3 . Bearing in mind that this is an upper estimate, and comparing 
to the bag model equation of state, e bag (p) = 4B + 3p, this result is still considerably larger 
than estimates with commonly assumed values of the bag constant B. 

Coming back to the region of the phase space where the quasiparticle model is well 
grounded, we finally address the question of the behavior of the pressure and the energy 



7 This value renders more precisely the rough estimate which was about 40% larger. Based on the 
pragmatic extension of the continuum extrapolation of the lattice data shown in Fig.n near T c , the fit 
led to a similar value for T s , but to A w 11. This demonstrates that details of the underlying lattice data 
are important for quantitative predictions at fi ^ but, on the other hand, that the estimates are rather 
robust. 
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FIG. 5: The estimate for the equation of state of quark matter at T = 0. 

density along the critical line near /i = 0. In the lattice simulations [^] both quantities have 
been found to be constant within the numerical errors. This is compatible with our result 
for small /i, 

p(T c (jj.),fi,) -p(T c (0),0) w -0.02 /i 2 T c 2 . (13) 

The corresponding change in the energy density is about three times larger. These results 
differ notably from the estimate from the bag model which, although the critical line has a 
similar shape for small /i, would yield coefficients larger by a factor of four. 



V. CONCLUSIONS 

Within our quasiparticle model we analyze recent n/ = 2 QCD lattice calculations 
of the equation of state at nonzero temperature and /i = 0, and then extend the quasiparticle 
model to nonzero baryon density. The resulting elliptic flow equation for the coupling relates 
the thermodynamic potential along the characteristic curves in the \i T plane. We argue that 
the characteristic line through T c {\i = 0) is related to the critical line in the phase diagram. 
This is confirmed by comparing our results for the curvature of the critical line atu = 0, and 
the variation of the equation of state along it, with recent lattice simulations [2| exploring 
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the region of small fi. 

We give an estimate for the equation of state of cold quark matter. Energy density and 
pressure are almost linearly related, as in the bag model, however with parameters obtained 
from the lattice data at \i — 0. The relevant physical scale is given by the transition 
temperature T c , and the parameter corresponding to the bag constant turns out to be large 
compared to conventional estimates, ~250MeV 4 . 

We have restricted ourselves to the case nj = 2, for which the lattice data for the equation 
of state at /i = appear to be established best. However, we expect similar results for other 
numbers of flavors since the pronounced decrease of the ratio p/po as T approaches T c , which 
indicates a large coupling strength, seems to be generic. This universality is then echoed at 
nonzero fi because for all nf the flow equation behaves for strong coupling similarly as in the 
perturbative limit, where a s is constant along the elliptic-like characteristics. Indeed, the 
shape of the phase boundary calculated in 19| for the physically relevant case rif = 2 + 1, 
although now being a crossover near T c , is very similar to the shape for n/ = 2. With 
the same reasoning, we remark that our estimates are robust with respect to remaining 
uncertainties of the underlying lattice data. Indeed, the equation of state at /i ^ is not 
very sensitive to the precise values of the model parameters as long as they reasonably 
describe the gross features of the equation of state at /z = 0. Therefore, the large energy 
density at smalpressure seems to be a general feature of the equation of state. 

As shown in this would allow for pure quark stars with masses < 1M Q and radii < 10 
km. Similar small and light quark stars have also been obtained within other approaches, 
cf. j^J. Such objects are of interest in the ongoing discussion of the data of the quark 



star candidate RXJ1856. 5-3754 



2l\ . It should be emphasized, however, that the outermost 



layers of such pure quark stars are metastable with respect to hadronic matter with a larger 
pressure at \i ~ n r . The details of the star structure depend sensitively on the hadronic 
equation of state j3| . However, as discussed in , a stable branch of hybrid stars with a 
dense quark core and a thin hadronic mantle could indeed be possible. 
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